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ABSTRACT 

>- ■ This is a preliminary report on the application of Difference Image Analysis (DIA) 

to galactic bulge images. The aim of this analysis is to increase the sensitivity to the 
! detection of gravitational microlensing. 

We discuss how the DIA technique simplifies the process of discovering microlensing 
events by detecting only objects which have variable flux. We illustrate how the DIA 
technique is not limited to detection of so called "pixel lensing" events, but can also 
be used to improve photometry for classical microlensing events by removing the effects 
Q-f of blending. We will present a method whereby DIA can be used to reveal the true 

Q ■ unblended colours, positions and light curves of microlensing events. 

We discuss the need for a technique to obtain the accurate microlensing time scales 



from blended sources, and present a possible solution to this problem using the existing 
HST colour magnitude diagrams of the galactic bulge and LMC. The use of such a 
solution with both classical and pixel microlensing searches is discussed. 

We show that one of the major causes of systematic noise in DIA is differential 
refraction. A technique for removing this systematic by effectively registering images 
to a common airmass is presented. Improvements to commonly used image differencing 
techniques are discussed. 
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1. INTRODUCTION 

The study of microlensing has become well es- 
tablished after a number of groups followed up the 
ground breaking proposal of Paczyhski (1986). Over 
the past few years well over a hundred events which 



have been discovered by the MACHO ( 


Alcock et al. 


1998a), EROS (Bcaulicu ct al. 1995), DUO (Alard 


k Guibcrt 1997 


) and OGLE ( 


Paczyhski et al. 1994) 



groups, in their attempt to characterise the nature of 
the dark matter halo of our galaxy. 

Fits to the light curves from microlensing events do 
not give us a direct measure of the mass of the objects 
causing the magnification. What they do provide is a 
measure of the lensing time scale and amplification. 
The amplification of light from a point source due 
to microlensing is given by equation |lj where u(t) 
can be obtained from equation ^, u m i n is the impact 
parameter in terms of Einstein radii and t max is the 
time at which maximum amplification is reached. 



A(t) = 



u (t) = u mm + 



u ( u 2+ 4)V2 

2(£ tmax} 



(i) 

(2) 



The time scale t of a microlensing event is set by 
the time the Einstein ring takes to traverse the source 
star at velocity V±. The projected size of the Einstein 
radius Re (equation ||), is dependent upon Dd, the 
observer-lens distance, D s , the source-observer dis- 
tance, Dd s , the lens source distance, and the lens mass 
M (in solar masses). 



Re 



4GM D d D 



ds 



(3) 



leads to increased accuracy in the determination of 
lens mass distribution. 

The number of events detected in microlensing sur- 
veys is limited by the number of stars which it is pos- 
sible to monitor and the overall detection efficiency 
In terms of telescopes, this means such surveys are 
limited by the size of the field-of-view of the survey 
telescope and its light gathering power. Limitations 
to such surveys also come from sampling rate and 
the seeing. These effects limit the accuracy of the 
photometry obtained and thus the event detection 
threshold. To maximise the number of events one 
naturally chooses to observe fields with the greatest 
density of stars. This introduces crowding which sets 
a seeing dependent limiting magnitude for such sur- 
veys. To maximise the number of monitored stars 
one can take an image in the best seeing conditions, 
where the crowding and sky background levels are at 
minimum. The detected stars can then be monitored 
even when the seeing is poor and sky background is 
high, by transforming the coordinates of stars from 
the reference image to each subsequent observation. 
This is the standard approach taken in such surveys. 

The use of crowded fields limits the accuracy of 
photometry. For images with poor seeing condi- 
tions, faint stars become immersed in the flux from 
neighbouring brighter stars. An important consider- 
ation when working with such fields, is that only a 
small fraction of the stars are actually detected. In 
fact, within each seeing disk there are generally many 
stars. The determined stellar centroid is actually a 
flux weighted mean centroid of the seeing disk. Any 
one of the stars within such a blend can be gravi- 
tationally lensed. Blending thus has a major effect 
on the actual number of stars that are monitored 
and hence the optical depth to microlensing. Such 
blending of microlensing events causes an amplifica- 



To uniquclv determine the mass of the lens caus- 


et al. 1993, Nemiroff 1994, 


Han 19971, 


Wozniak & 


t anv given event, one must determine its distance. 


Paczyhski 1997, Alard 1997 


. In past analyses such 



the source distance and its transverse velocities rela- 



effects have been taken into account statistically (Al 



tiv a to our linc-of-sight. This is usually not possible. 



However, for a small number of exotic events where 
binarity or parallax are found, the degeneracy of M, 



cock et al. 1997a ), but correction s are difficu lt to 
apply on an event by event basis ( Han 1997b ). In 



Dd/D s and VI. (e.g. Alcock et al. 1995) can be bro- 



ken. One can instead extract the distribution of lens 
masses by assuming the distribution of lens transverse 
velocities and distances ( |Griest 1991 ) . As this is a sta- 
tistical process, an increase in the number of events 



this way the method of monitoring a star field using 
fixed positions, in some sense, does allow us to mon- 
itor some sources which were too faint to have been 
detected by virtue of their own flux. This does not, 
however, allow us to detect any microlensing events 
due to such source stars which are not blended with 
a neighbouring star or stars bright enough to be de- 
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tected. 

New, large, ground based telescopes are useful in 
increasing the number and frequency with which star 
fields can be monitored because of shorter exposure 
times are required. But such telescopes are still prone 
to the same intrinsic seeing limits as small telescopes, 
and thus the same crowding limits of images. The ef- 
fec tive seeing can be improved with the use of better 



obsE 



Tving sites, but this does not resolve the issue of 



blending completely when many stars are present per 
square arc second. Attention must be applied to im- 
proving the reduction techniques to remove blending 
effects. Such techniques would also ideally increase 
the number of events detected by increasing the effec- 
tive number of stars monitored per image. To achieve 
both requires a technique which will overcome the 
crowding in the images and increase the monitored 
area within each image. With this intent we shall ap- 
ply difference image analysis (hereafter DIA0). This 
uses a similar technique to that first described by 
Crotts (1992) and Phillips & Davis (1995). 

Lensing events where the source stars are only de- 
tectable during microlensing are usually designated 
'pixel lensing' events. Events where the source star is 
a re solved are term ed 'classical micr olensing' events 
(see pould (1996)1 and pould (1997)] for full details). 
The term pixel lensing comes from an analogy with 
microlensing surveys in M31 (Ansari et al. 1997), 



where each pixel in this survey represents the sum of 
hundreds of unresolved stars that may be microlensed. 
This term is applicable with local microlensing sur- 
veys (bulge, LMC) to an extent, since even for the 
bulge, given arc second seeing, each pixel contains 
light from a number of stars. In this way we define lo- 
cal pixel lensing events as microlensing events where 
the sources are not associated with stars we moni- 
tor. Furthermore, to remove the ambiguity in what 
is a pixel lensing event and what is a classical lensing 
event, we will use the colour and position information 
available with the DIA technique. The three proper- 
ties flux, colour, position allow us to make a fairly 
robust separation between these event types. 

In the following sections we will demonstrate how 
the DIA technique can be used to increase the num- 
ber and quality of results from microlensing surveys. 
In the next section we will briefly outline the observa- 
tional strategy. In §3 we will also discuss how the DIA 
technique was employed on our set of data and what 



improvements to the standard method were made. In 
§4 we will discuss the relevance of DIA to the issue of 
blending. In §5 a comparison between PSF photome- 
try and this technique is made. In §6 we outline how 
to determine source fluxes from the results. In the 
final section we shall make our concluding remarks. 
The light curves and parameters for the microlensing 



events discovered will be presented in ALcock et al 



(1999 paper2) 



2. OBSERVATIONS 

The MACHO observation database consists of over 
70,000 individual observations of the galactic bulge, 
the LMC and the SMC. The present analysis con- 
siders only a single 42' by 42' field (Macho id 108) 
in the galactic bulge centred at a = 18°01'20", 
S = -28°17'39" (J2000). Observations were taken on 
the Mount Stromlo and Siding Spring Observatories' 
1.3m Great Melbourne Telescope with the dual-colour 
(red, blue) wide-field "Macho camera" . All bulge ob- 
servations have 150 second exposure times. 

The Macho camera consists of a mosaic of eight 
2k x 2k CCDs (four red, four blue) with 0.63" pixels. 
Each of the eight CCDs in the Macho camera consists 
of two amps. Light passing through the telescope is 
separated into two passbands using a dichroic beam 
splitter allowing both red and blue images of a field 



'This technique is also referred to as image subtraction. 



to be taken simultaneously (Marshall et al. 1993). 

Observations in this analysis span the dates from 
the 10th of March 1995 to the 2nd of August 1997 
with breaks from the end of October to the beginning 
of March each year when the bulge is unobservable 
at Mount Stromlo. Within this observing period 385 
observations of the target field were taken. The seeing 
for the data set varies from 1.2" to 6.5" with a mean 
of 2.3". The sky background level in the blue band- 
pass varies from 1000 to 30000 counts with a mean 
~ 2600 counts. Similar levels were observed for red 
images. 

We chose to reject a number of observations from 
this dataset as good photometry was required to de- 
tect of microlensing of faint stars. Images with seeing 
FWHM > 4" or blue-band sky background > 8000 
ADU were excluded from the reduction. These two 
cuts removed 42 observations from the data set. As 
the usable area of a difference images is dependent 
on the differences in pointing between observations, 
we also excluded a small number of images where the 
difference in pointing of the reference and subsequent 
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observation was greater than 25". With this criteria 
imposed we rejected 19 more observations leaving 324 
before reduction was attempted. 

3. DIFFERENCE IMAGE ANALYSIS 

Since the work of Crotts (1992) and Phillips & 
Davis (1995) there have been a number of applica- 
tions of DIA type techniques (Tomaney & Crotts 



1996| , [Rciss et al. 1998| , |Alard fc Lupton 1997|) . These 



papers show the rapid evolution of the technique to a 
near optimal case. The implementation of the tech- 
nique as used here varies in several important aspects 
to these approaches. So we shall outline our tech- 
nique, noting the similarities and differences to these 
applications. 

In brief, DIA of our images involves registering a 
test image to a preselected, high signal-to-noise ratio 
(hereafter S/N), low airmass, stacked reference image. 
Next ~ 200 bright, uncrowded, so called PSF stars, 
are selected to calculate the convolution required to 
map the template image to the test one. The convo- 
lution kernel for this mapping is calculated and ap- 
plied to the reference. The test image is then pho- 
tometrically normalised (scaled and offset) to match 
the convolved reference image. The resultant, regis- 
tered, convolved and normalised image is then sub- 
tracted from the test image. The process is carried 
out for red and blue passbands, and objects are de- 
tected. The positions of these objects are matched in 
the two colours to separate real objects from spurious 
ones. The results are then sorted and characterised 
to separate microlensing events from variable stars. 

3.1. Template construction 

To perform DIA the first requirement is to have a 
comparison reference image (template) to difference 
against. To minimise the noise contribution to the 
final difference images from such a reference image, it 
is an obvious step to use the highest S/N image. The 
MACHO project has a database of hundreds of obser- 
vations of each field taken with the same telescope, 
under similar observing conditions. This database 
makes it possible to increase the S/N for a reference 
image simply by stacking matched images. As the 
the highest S/N is associated with the best seeing, 
we choose to degrade the reference image seeing to 
the test image. Since we require that the reference 
image has good seeing we can only combine images 
with good seeing to create it. This considerably lim- 



its the number of images we can combine to make the 
reference image. 

To produce this reference image we require the con- 
stituent images be taken at low airmass. This is re- 
quired because the combination must be unaffected 
by refraction (see section 3.3). We also require that 
such images do not have high sky backgrounds or sky 
gradients, so that the combined sky level contributes 
little noise to the final difference images. Low sky 
levels allow us to place the lowest possible detection 
thresholds for pixel lensing events. Care has to be 
taken selecting images, as even a small residual sky 
gradient in the combined reference image can affect 
how well two images are photometrically normalised 
and thus the quality of the entire set of difference im- 
ages produced. 

Aside from the advantage of an increased S/N in 
the combined reference image, the process of stack- 
ing images has the effect of removing most of the bad 
pixel regions in the images. Small features in the 
images due to differences in sky background, trans- 
parency, fringing, or poor flatfielding are also reduced. 

The airmass, seeing, and sky level selection criteria 
strongly restrict the number of observations suitable 
for the reference image combination. Among our ~ 
324 observations, only five meet all our requirements. 

3.2. Image Registration 

The importance of obtaining an accurate registra- 
tion when performing DIA can not be over empha- 
sised. To accomplish an accurate registration, and 
to reduce the time required for analysis, we took ini- 
tial estimates for the geometric alignment transfor- 
mations from a database produced when each image 
was PSF photomctcrcd. These transformations were 
applied to some hundreds of bright stars spread across 
each image. The transformed positions of these stars 
were then used as starting points where more accu- 
rate centroid positions were determined. The initial 
centroid positions were found to be good to within a 
pixel (0.63"). Such accuracy means we do not have 
little possibility of centering on the wrong star. These 
new centroid positions for the stars in the test and 
reference images were used in the IRAF tasks geomap 
and geotran to register the image. The use of a large 
number of stars allowed us to determine an accurate 
mapping between the images where small rotations, 
translations and distortions may be present. 

The images were initial registered with a simple 
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linear transformation on a 10' x 10' scale. However, 
we found this was inadequate for the alignment ac- 
curacy we require from this technique (< 0.1"). This 
is likely to be caused by distortions between images, 



due to small effects such as differences in telescope 



flexure with zenith distance. For the final analysis we 
have used a fifth order Chebyshev geometric trans- 
formation mapping for blue images, and a third order 
mapping for red images (where differential refraction 
effects are not so important, see next section). 

The pointing accuracy of the Great Melbourne 
Telescope is good, with the standard deviation in off- 
set between successive observations of a field being 
approximately 10". For difference images this results 
in the loss of a very small amount of data at the edges 
of the images where the observation and reference im- 
age do not overlap. 

The registration process is accurate to abou t 0.07 
pix els (~ 4 5 milli-arc seconds) on average (c.f. Fteiss 
et i d. 1998] , ~ 0.3"). We believe our accuracy is due 



to the inclusion of higher order transformation terms 
and the large number of stars used to constrain the fit. 
To show the importance of accurate image registra- 
tion we present figure [l]. In this figure we approximate 
a standard stellar profile as a Gaussian. Using the 
average seeing of ~ 2.5" and an offset between two 
stellar profiles of 0.3", we expect the average resid- 
ual (as seen perpendicular to the offset direction) to 
have the form of the long dashed line given. For our 
average registration accuracy (0.07") the stellar and 
residual profiles are given by the short dashed lines. 
For this figure we did not include the effects of pixeli- 
sation or photon noise, as this makes little difference 
to the observed effect. 

3.3. Differential Refraction Corrections 

The first mention of the deleterious effects of atmo- 
spheric refraction on difference images was made by 
Tomaney and Crotts (1996). They stressed that, for 
broad-band filters, the centroid position of a star is 
dependent upon the airmass of the observation. The 
phenomenon is well known in astrometry and has the 
effect of increasingly elongating the PSF with air- 
mass. For observations taken at different airmasses 
this leads to poor subtraction. By poor subtraction 
we mean that the systematic noise contribution from 
the refraction effect is obvious in the images. The 
presence of this is revealed by a significant residual 
flux remaining in the difference image. 



To better understand the refraction effect, one 
must understand that every star within the frame has 
a centroid position which is dependent upon its colour 
and the airmass at which it is observed (Gubler & 
Tytler 1998). If all the stars were of the same colour 



then the registration process alone would compensate 
for the shift in centroid position with airmass, and 
the difference image would only be effected by the 
elongation. However, this is not the case. The image 
registration is in reality a registration to the average 
colour of the stars. We suggest that, two images ob- 
served through broad-band filters at blue wavelengths 
and at moderately different airmasses, can not be ge- 
ometrically registered accurately enough to provide 
truly photon noise limited difference images. How- 
ever, we further suggest that, if each star were to 
be registered separately we could come close to this. 
With this realisation in mind, Tomaney and Crotts 
(1996) proposed that 'There was no easy way to com- 
pensate for this, particularly without knowledge of 
the precise position and colour of all detected and 
undetected (blended) stars'. Similar effects have also 
been noted by Melchior et al. (1996) who believed 
that atmospheric dispersion affected their pixel lens- 
ing search in M31. This effect was given to be be- 
tween 0.73" to 2.75" in the extremes of their blue 
filter. With this result it is easy to understand the 
importance such an effect could have. However, such 
a large effect does not occur in our DIA analysis, even 
in the worst the first order refraction effect 

is an offset of all stars in an observation. This offset 
is taken out in the registration process. It is only the 
smaller second order effects, caused by differences in 
star colour which cause differential refraction offsets. 
The corrections for these effects are in the order of 
tenths of an arc second at blue band wavelengths and 
so are still appreciable. 

To correct for differential refraction the following 
approach has been taken. As a first approximation, 
the effective temperature for any given star can be 
associated with a colour. We know the position of the 
centroid of a star is dependent upon the colour of the 
star and the airmass. So if one knows the airmass, the 
colour and the observation parallactic angle it should 
be possible to determine the change to the stellar flux 
distribution with airmass and parallactic angle. 

To determine the degree of change we first mod- 
elled our galactic bulge stars based on Bessell's ( 1998 ) 
tabulated colour-temperature relations. We then as- 
sumed a blackbody spectrum approximation for these 
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stars with no blanketing or spectral features. 

These spectra were then integrated within the MA- 
CHO Camera's blue(B m ) and red(i? m ) passband re- 
sponses to find the average centroid wavelength asso- 
ciated with a given star colour. The V-R colours for 
Bessell's (1998) tabulated stars were transformed to 
Macho R m and B m colours using the results of Alves 
(1998). A star of a given colour can thus be related 
to this centroid wavelength. 



Rn — 



R 



in 2 - 1) 
2n 2 



Zt ~ z a 



R » R Q tanz a - 0".067tan 3 z 



(4) 



(5) 



(6) 



The refraction of light at a given wavelength, tem- 
perature and airmass is given by |Filippenko (1982)| . 
Using the equations in Filippenko (1982) and equa- 
tion ^ one may determine the constant of refraction, 
where n is the refractive index of air. The degree of 
refraction is given by equation^, where z t is the true 
zenith position and z a the apparent. Under normal 
temperature and pressure conditions (15°C, 760 mm 
Hg) this reduces to equation |(| and Rq becomes 58.3". 

The positional shift of the centroid of the star is 
associated its average blackbody wavelength. For a 
star of known colour the centroid shift can be ob- 
tained. The offset in the centroid positions for the 
model stars are given in figure ^ for three values of 
airmass. To make use of this information, one needs 
to know the colour of each star and must be able to 
shift its centroid in proportion to this colour. This 
task sounds more daunting than it really is. In fact, 
this can easily be achieved (at least approximately), 
for two images taken at a known airmass in the follow- 
ing way. Using the IRAF tasks geomap and geotran, 
it is possible to map two images of different sizes, 
orientations and even geometric distortions onto each 
other. We used these tasks to map the images from 
the two passbands of the MACHO camera onto each 
other. These mapped images were then used to form a 
quotient image which, after sky background subtrac- 
tion, can serve as a colour map of each pixel in the 
field. However, as the centroid wavelengths of each 
star within the two passbands varies with airmass, 
one must calculate this quotient near unit airmass not 
to have a colour map which is affected by refraction. 



To find the relationship between the pixel values 
of our colour map, and the V-R colours of the stars 
present in the colour map, a simple calibration is car- 
ried out. Having performed photometry on the indi- 
vidual images in the colour map, one can associate a 
B m — R m and hence, a V-R colour, with a quotient 
image pixel value. We thus have a V-R colour for 
every pixel within the image (except for saturated or 
bad pixels). The relationship between the measured 
B m — R m photometry values and the associated colour 
map pixel value is given in figure [|. The relationship 
is quite strong for most stars. A small number of 
points are scattered due to stars being blended with 
neighbours. 

As the degree of refraction is wavelength depen- 
dent it is in fact only necessary for us to apply our 
corrections to images taken through the blue pass- 
band. In red passband the centroid offset is less than 
one quarter (< 0.1") that of the blue passband (see 
figure |. 

The final result is that each blue image pixel value 
in the reference image was interpolated based on the 
calibrated colour map pixel values, the airmass of the 
observation relative to the reference image, and the 
parallactic angle of the observation. The flux correc- 
tions were carried out while imposing a condition of 
flux conservation within the image. We interpolated 
the flux in the reference image rather than in each 
observation as, for most observation it is not possible 
to make a colour map for these due to first order of 
refraction effects. Colour maps made with individual 
observations would have much lower S/N and CCD 
defects. 

The results from this reduction were used as an 
initial estimate, as it was known that the results for 
the differential refraction offsets would be dependent 
on the assumptions about the model stars (smooth 
blackbody emission) and the uncertainty in the form 
of the Macho blue band response function. This es- 
timate was improved by differencing a number of im- 
ages taken at a range of airmasses to obtain a semi- 
empirical result for the offset required with colour and 
airmass (relative to the reference image). 

The scale of our corrections can be seen in figure 
§. On the left is a difference image where the effects 
of differential refraction have not been removed. On 
the right the differential refraction corrected images 
is shown. The test image used in producing these dif- 
ference images is at an airmass of 2.4 (approximately 
the airmass limit imposed on these observations) . The 
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template airmass is 1.01, so the corrections applied in 
this case are quite large. However, the same degree of 
differential refraction effects is observed at much lower 
airmasses when the seeing quality of the test image 
is good. The residual is highly dependent on seeing, 
with the best seeing images being affected much more 
because of the steepness of the stellar profile. 

The calibrated differential refraction corrections 
for the blue passband of the Macho camera, as a func- 
tion of airmass and colour, are given by: 



S = 0.29 x [tan(acos{l/A R ) - tan(acos(l/A T ))} (7) 



= Sx (C-0.6) 



(8) 



S is related to the scale of the offset for airmasses 
An and At , of the reference and test images respec- 
tively. O is the offset to apply to a given pixel (in arc 
seconds) of V-R colour C. This process compensate 
for the shift in the stars centroid position. It does 
not compensate for the dispersion of the PSF with 
airmass. However, we believe this effect to be quite 
small and for our purposes ignorable. 

3.4. PSF Matching 

The Point Spread Function (PSF) shape of two im- 
ages taken at different times is never exactly the same. 
If we were to match the photometric conditions of a 
reference and test image and then difference the im- 
ages, we would not expect to obtain a resulting image 
without systematic noise. What one would invariably 
find from such a process is that at the position of each 
star there would be systematic residuals. The degree 
and structure of these residuals would be dependent 
upon the form of the spatial difference between the 
PSFs. To achieve the best difference images one has 
to match the form of the PSFs. As we mentioned ear- 
lier, the template is constructed so that it has high 
S/N and seeing similar to the best observation^. To 
match the PSFs we degrade the reference images' see- 
ing to match that of the test observation. 

The PSF-matching process is based on the fact 
that it is theoretically possible to match the profiles of 
stars observed under two different sets of conditions, 
wit h a oimplc convolution of the form given by oqua 



tioi l (Q). — Hero r charactorioca the flux diotribution 



of a star in good (better) seeing and t that in poor 
seeing, k is the convolution kernel. A convolution in 
real space is equivalent to a multiplication in Fourier 
space. Therefore, in principle, the Fourier Transform 
of the kernel required to match the good seeing (ref- 
erence) image to the poor seeing (test) image, should 
be the quotient of the Fourier Transforms of the star. 
The Inverse Fourier Transform recovers the required 
matching kernel (equation FLO) . 



t{x, y, z) = r(x, y, z) * k(x, y, z) 



\FT(r) 



(9) 



(10) 



In reality a division of these Fourier Transforms 
is very sensitive to the high frequency, low power, 
noise component. This leads to a poor match of the 
images. To fix this one can use the fact that, the 
PSF of a star is roughly Gaussian, and the Fourier 
Transformation of a Gaussian is a Gaussian. We can 
then select a level in Fourier space below which the 
noise component is dominant, and replace this with a 
Gaussian fit to the FT (Ciardullo 1990). This method 
is not always useful because in many cases the wings 
of the FT are not well modelled by a Gaussian. 

To determine the best convolution kernel, the high- 
est S/N, least blended PSF stars are required. We 
thus need these PSF stars to be unblended with neigh- 
bouring stars in the bulge where effectively all stars 
are blended. Blended stars must be removed from 
those used in determining the PSF. 

This situation was overcome by producing a list 
of a couple of hundred bright stars on a given im- 
age. Stars were culled from the list based on the rel- 
ative proximity and brightness of neighbouring stars. 
Blending of bright stars can in part be determined 
from the shape of the PSF of each star. We thus se- 
lected stars based on their ellipticity, position angle, 
FWHM and moment. The remaining stars were then 
combined to form a generalised PSF profile for the 
image. 

A further complication to the ideal case arises from 
the fact that the form of the PSF varies across a frame 
because of poor focussing, telescope flexure and tem- 
perature dependent effects (see Tomaney & Crotts 
1996| ). If a single star is used to PSF-match an entire 



15 A11 observations are well sampled with minimum FWHM of 
~ 2.8 pixels. 



image, the systematic noise in the difference image in- 
creases the further one gets from that star. A solution 
to this is to split the image into sub-rasters and use 
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a local PSF to match the sub-regions of the image. 
The sub-rasters can then be mosaiced back together 
to reform the image. Thankfully the magnitude of 
the PSF variation for our observations is small, with 
modest sized effects due to differences only being seen 
on the scale of 500 x 500 pixels (or 5' x 5'). 

The actual matching of PSFs is accomplished us- 
ing the IRAF task PSFmatch. This task is similar to 
the implementation used by Riess et al. (1998) but 
has been updated to include new features. The con- 
volution kernel is determined for sub-rasters 500 x 
500 and each is convolved separately then mosaiced 
to form lk x lk images. The entire observation field 
of the MACHO camera is 4k x 4k per bandpass but 
larger images were not made because of problems with 
matching such as the scaling, colour terms, di fferences 
in bias le vels and gain in the different CCDs ( Reiss et 



al. 1998) 



The process of replacing the noise with the Gaus- 
sian fit is not regarded as the optimal approach be- 
cause of there are real differences between the PSF 
shape and a Gaussian. To improve the difference im- 
age quality we characterise the residuals in the differ- 
ence image and remove them. This is accomplished 
in the following manner. By stacking sub-rasters of 
the difference image at the positions where the PSF 
stars were in the reference image, we determined the 
median systematic residual at each pixel of the PSF 
profile. This stacked residual image characterises the 
real deviations of the PSF shape from that used in the 
matching process. This image was then scaled and 
subtracted from our test image PSF. We finally had 
an essentially empirical PSF and no need to replace 
the wings of the FT. The convolution kernel is now 
recalculated without the Gaussian replacement and 
the difference images recalculated. This final PSF ap- 
peared to give the best possible subtraction with the 
available stellar information. This step was found to 
reduce the systematic residuals by a factor of ~ 2 in 
most cases. The final resulting difference images had 
an average systematic noise component of ~ 1.3% in 
both colours, with some images having less than 0.5%. 

3.5. Photometric Scaling 

A standard approach to matching reference and 
test images photometrically is to make a linear fit to 
the pixel values for a single star found on both images 
(Reiss et al. 1998). The test image is then scaled 
and offset using the fit values. We refer to this as a 
single point calibration. Another approach is to use 



photometry taken on a matched set of stars in the two 
images and scale the image based on this information. 
In our analysis we apply this second approach which 
we believe is superior. 

The single point technique has the advantage that 
it uses all the pixel information for this region. How- 
ever, it suffers from the fact that most of the pixels 
have a low S/N ratio as they come from the wings 
or background around a star. As there are very few 
pixels with high S/N, the fitted slope can be skewed 
by the pixels near the background noise limit. 

An important consideration when matching images 
is that the two images were probably taken under 
different seeing conditions. Although the PSF pro- 
files have already been matched, real images can have 
gradients in sky brightness and transmission. Sky 
brightness gradients can often be attributed to the 
proximity of the observed field to the moon, and nat- 
urally are dependent on its phase. Differences in the 
transmission come with airmass, cloud cover and dust 
extinction. Such differences in sensitivity can also be 
caused by problems with the flatfielding. 

A single point calibration can not compensate for 
the presence of these spatial variations. However, this 
can easily be achieved by performing photometry on 
matched pairs of stars across the two images. This 
photometry gives a scaling factor (which represents 
changes in transmission), and an offset term (which 
represent the change in the sky brightness), for each 
position where the photometry was performed. To 
deal with the possibility of gradients we fit a low or- 
der polynomial^ to the scaling and offset terms and 
determine a transformation for each pixel in the test 
image. Such a method requires a large number of 
points within the image to constrain the fit. This is 
not a problem as there are many thousands of stars 
per square arc minute in bulge fields. 

The photometric scaling we used was always ap- 
plied to the test image so that images were always 
registered to the same reference image. In this way 
no correction was required for transmission differences 
with airmass. 

3.6. Object Detection 

The fully matched images were differenced to re- 
veal objects whose flux rate has varied in some way. 
The images were then searched for these sources. 



16 We expect variations in sky level and transmission to be rela- 
tively smooth. 
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Variable stars are detected as positive or negative 
sources as they became brighter and fainter than they 
were in the reference image. However, microlensing 
events and asteroids generally only appear as positive 
excursions from the reference image baseline^. 

For source detection we applied a purpose written 
programme. We felt this approach was necessary as 
the nature of the noise distribution is unlike that in 
other images. The fundamental difference is that the 
noise varies from pixel to pixel even though there is 
no signal. This is due to the photon noise of the 
subtracted stars. With such a programme we search 
for positive and negative source simultaneously. 

To detect objects we would ideally would like to 
know the noise at each point in the image. To de- 
termine this we simultaneously examine the reference 
and difference images. The significance of each detec- 
tion is accessed based on its S/N. We calculate the 
noise at each source position using the flux in the 
reference image and systematics from the difference 
image. The various noise components (photon, read- 
out, systematic) are well characterised in the reduc- 
tion and used to provide a position dependent noise 
threshold for the difference image. The photon noise 
of the test image dominates for faint objects whereas 
for brighter objects the systematic noise from the dif- 
ference image becomes important. As the seeing in 
our data set is always greater than two pixels, we 
choose to determine the signal in a 3 x 3 pixel box 
at every point in the image. The amount of signal in 
this box is compared with the noise expected for the 
same pixels in the reference image. Results with a 
S/N > 3 were written to a file along with other noise 
parameters such as the systematic and readout noise 
of the difference image. 

The detection process is carried out with images 
from both passbands. The positions of objects found 
in the two colours are taken to be matched if they 
were within 1". This colour match provides a strong 
constraint which separates the real objects from the 
large number of detections due to cosmic rays, bad 
pixels and saturation effects. 

Candidate objects which pass these selection crite- 
ria were then checked against those obtained in previ- 
ous reductions of the field. If no match is found, these 
new results were added to the database of object po- 
sitions. If a result has previously been detected, the 



17 It is possible for microlensing events to appear with negative 
flux if they were amplified in the combined reference images. 



number of detections for this object is incremented. 
This provides us with an extra parameter in char- 
acterising the nature of the object. Objects detected 
most often can usually be attributed to variable stars, 
whereas other detections can be due to asteroids or 
microlensing events. 

3.7. Light Curves 

Aperture photometry is performed on the red and 
blue difference images at the sources positions. We 
feel PSF photometry is unnecessary as the difference 
images are no longer crowded with non-variable ob- 
jects. The noise is obtained from the initial images as 
noted in the detection strategy. With this photome- 
try we produce a database which is used to produce 
difference-flux light curves for microlensing event de- 
tection. 

These positions where photometry is performed are 
independent of whether there is flux at this position in 
the reference image or not. This varies from the tradi- 
tional approach which has been to only follow up stars 
detected in the template (reference) image. With this 
technique we are searching among a database of vari- 
ables for microlensing rather than among a database 
of all objects. The number of light curves which re- 
quire scrutiny is thus around a factor of twenty less 
than the number of stars detected in the reference 
image. The removal of the positional dependence 
gives us greater sensitivity to detecting microlensing 
of stars which were either too faint or too blended to 
be resolved in the reference image. This effectively 
provides us with a larger search area and a greater 
number of monitored stars. 

4. BLENDING 

Firstly, let us define what we mean by blending. 
This is the case in which two or more stars are lo- 
cated within a few pixels of each other and the see- 
ing disk is greater than a few pixels. Such blending 
is dependent on seeing and masks the true baseline 
flux of microlensing sources. In this analysis, we are 
not concerned with the case where the lens in a mi- 
crolensing event emits flux which is blended with the 
flux from the source. In this case the flux of the lens is 
not amplified and such blending is generally negligi- 
ble for stars with masses consistent with microlensing 
results. The colour shifts associated with this type of 
blending usually can not be detected with the present 
level of photometry (Buchalter et al. 1996). This case 
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is also generally inseparable from blending with ob- 
jects which are not the sources. We are also uncon- 
cerned with blending where the source is in a binary 
association. Here the sources can be blended within 
the Einstein radius, but this constitutes a small num- 



ber of events (Griest fc Hu 1992, Dominik 1998) and 



in many cases one can detect such events from the 
microlensing light curve. 

The presence of the apparent blending of the stars 
monitored in microlensing surveys has been known 
for some time (Nemiroff 1994), and is inherent in the 
the fact that the fields used for microlensing searches 
are crowded. These fields were chosen because of the 
large number of stars that could be monitored at one 
time. Crowding is demonstrated by figure |^, where 
the number of stars blended within the seeing disk of 
the bright stars can be clearly seen. 

The main effect of the blending is that we do not 
know the true source flux for most events. If such 
events are fitted as unblended sources, the ampli- 
tudes of the events and the event time scales are un- 
derestimated, and the number of stars monitored is 



under-estimated ( Han 1997a , Wozniak fc Paczyhski 
19S0). The overall consequence, if blending is unac- 
counted for, is that ~ 40% of events towards the bulge 
are affected by amplification bias and thus the deter- 
mined optical depth can be overestimated by a factor 
of w 1.3 - 2.4 ( |Han 1997a| ). The blending effect is 
therefore, a major factor in the determination of the 
true optical depth to microlensing toward the galactic 
bulge, LMC, SMC and M31. Statistical corrections 
have been made for the LMC and the bulge optical 



depths ( [Alcock ct al. 19974 |Alcock ct al. 1997b| ), but 
corrections are really required on an event by event 
basis. One way to bypass the effects of blending would 
be to have HST images of all the events. With this 
approach we could resolve the source star and obtain 
an unblended source flux. Unfortunately a large num- 
ber of hours of HST time are needed to perform this 
( Han 1997b ). We note, however, we do have HST 
observations of microlensing sources for a number of 
important events. 

To determine the blending in the case of classi- 
cal microlensing, one must attempt to fit the light 
curves to find the unlensed, blended, flux component 
of the source baseline flux. In this regard Wozniak & 
Paczyhski (1997) noted that one can only determine 
the blend fraction, with the present accuracy of pho- 
tometry, when the impact parameter is small, u m i n < 
0.3. This is the case in only a small percentage of 



events. The situation is not quite the same with DIA 
as we have different information. We shall outline this 
in the next section. In short, instead of the colour, 
centroid and light curve of the blend, we have the 
colour and centroid of the source and the light curve 
of the source amplification. 

5. DIA VERSUS PSF PHOTOMETRY 

Our usual technique to detect microlensing is to 
perform PSF photometry simultaneously at fixed po- 
sitions in two passbands. With our DIA technique, 
photometry is carried out at the end of the reduction 
process, at positions where excess flux was detected 
at some time during the reduction. 

The main advantage of using traditional PSF pho- 
tometry over DIA comes when dealing with fields 
where the stellar profiles are not blended together. In 
such situations the determination of the PSF from a 
few stars is relatively simple. PSF photometry, unlike 
DIA photometry, is not subject to either the addition 
of photon noise via the differencing process or the 
systematics introduced in the image alignment and 
matching processes. 

One common belief about difference image tech- 
niques, is that it is not as powerful at constraining 
microlensing as standard PSF photometry because it 
does not give a baseline flux. This is not the case. For 
classical microlensing events, photometry can simply 
be done once on the reference image to determine the 
baseline flux. In our present method of analysis this 
image has a higher S /N than any individual image in 
the data set, and has the best seeing and a well de- 
fined PSF. The photometry baseline determined from 
this frame can simply be added to the individual dif- 
ference frame photometry to provide the same flux 
baseline as the standard PSF photometry. However, 
this source flux baseline can still be blended as in the 
PSF case. 

To produce difference images with the lowest pos- 
sible systematic noise, the input images must be ac- 
curately photometrically matched. If the difference 
image has only a small contribution from systemat- 
ics, then we know that the images were accurately 
matched photometrically. To perform PSF photome- 
try we have to make corrections for airmass and dif- 
ferences in transmission between observations. 

In figure [?] we show the photometry performed on 
an object with DIA and with PSF photometry us- 
ing the same set of images. A quite dramatic im- 
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provement in the photometry is seen. Proof that this 
demonstrates a real improvement in photometry is ev- 
ident from the microlensing fit residuals shown on the 
same figure. Although only one colour is shown, data 
points from both passbands were consistent within 
the uncertainties presented. The improvement in pho- 
tometry comes partly from the fact that the source 
in this event is highly blended with a neighbouring 
bright star. With DIA the nearby non-variable star 
is more accurately removed than with the standard 
PSF photometry. The PSF photometry is also de- 
pendent on how well the centroid positions of the 
stars were initially determined in the template im- 
age. On the other hand, DIA is dependent on how 
well the entire frame was aligned. The flux from un- 
lensed blended stars is removed accurately even when 
unresolved with DIA. Such an improvement in the 
baselines of microlensing events are required for the 
detection of the subtle effects of parallax and the pres- 
ence of binary sources. 

Another advantage of the DIA scheme over PSF 
photometry is that the difference flux represents the 
true colour of a microlensing event. This colour in- 
formation is useful in deciding whether an event is 
blended or not. Furthermore, all difference images are 
matched spatially and photometrically. Thus we can 
stack images over the period where the microlensed 
flux was above the background noise. The stacked 
difference image can then be used to determine the 
source star colour more accurately. 

The present solution to the problem of blending in 
the PSF photometry is to use the fact that the mea- 
sured flux within a microlensing event can be repre- 
sented by equations [ll] and [l2] ( ALcock et al. 1997a), 



fR(t) = fuR + A(t)f, 



f B (t)= fuB+A(t)f sB 



(11) 



(12) 



where f U R and f U B are the fluxes of the unlensed 
blended sources and / s _r, J s b are the baseline fluxes 
of the lensed star. All of these terms must be found 
from the microlensing fits. The major problem with 
this approach is that, even if we know the true colour 
of the source, we do not know its true brightness. For 
bright, well covered microlensing events, source fluxes 
can be found. However, when uncertainties of > 1% 
are associated with the photometry, the fit is prac- 
tically degenerate with respect to the unlensed flux 
component (see Wozniak & Paczyhski 1997). This 



means we can not accurately determine the value of 
the amplification for a large number of events. 

In figure || we show three models of the amplitude 
of a microlensing event with different values of the 
lensing parameters corresponding to different values 
of source flux (given in table [l]). In the top panel there 
are three curves which are almost indistinguishable. 
In the lower panel we plot the difference between the 
dashed curves and the solid one. One can see that 
even for such large differences in source flux, hence 
lensing amplitude, differences of less than a couple of 
percent occur in the form of the difference in curve 
shape. The difference between the curves increases 
only slowly as one moves away from the true values. 
This predicament is improved greatly with photome- 
try taken in two or more colours because of the extra 
leverage this gives. 

For DIA the flux in the two band-passes is given 
simply by equations 13 and|IJ. Here we do not have 
blending to consider so there is no unlensed flux term. 
The difference fluxes in the two passbands are given 
by fnD(t) and fsDit)- The constant C is the ratio 
of blue to red flux which can be obtained from the 
colour of the source. 



f RD {t) = (A(t) - l)f sR 



f BD (t) = (A(t) - l)Cf, 



sR 



(13) 



(14) 



For this situation we only need to be able to mea- 
sure the baseline flux in one colour to determine the 
amplitude, since Cf s ji = f S B- If this is not possi- 
ble, the accuracy of the determination of the ampli- 
tude suffers from the same problems as the PSF pho- 
tometry. For cases where we can not determine the 
baseline flux (pixel lensing), it should be possible to 
estimate this quantity statistically. 

6. SOURCE FLUXES 

The case of determining the source flux with pixel 
microlensing is similar to the case of determining the 
blending for classical microlensing. In both situations 
we do not know the initial source flux. The standard 
way to determine the true amplitude of microlensing 
events is to use the shape of the light curve. The ac- 
curate determination of the source flux requires well 
sampled light curves with small percentage errors. In 
many cases the present data does not meet this re- 
quirement. 
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With pixel lensing the source star of the event is 
initially unresolved and the S/N is generally low. The 
task of determining an accurate source flux can be im- 
possible in many cases. However, from the reference 
image we can determine an upper limit to this flux. In 
such cases, as for DIA of classical microlensing events, 
we still have an accurate position and colour for the 
source. With this colour (and the associated uncer- 
tainties), the HST luminosity function of the bulge 
(Holtzman ct al. 1998), the upper limit of source 



flux, and the shape of the light curve, one can deter- 
mine the probability distribution associated with the 
source flux for each pixel lensing event. This source 
flux probability distribution and the light curve shape 
can be used to give the distribution of i for each event. 
The combination of these t distributions gives us an 
overall distribution of t from which the lens mass dis- 
tribution can be extracted in the traditional way (Gri- 
est 1991). 

For classical microlensing events where the sources 
are faint and the amplifications are moderate, the as- 
sociated photometry errors are typically greater than 
a few percent. Typically for these events, as for pixel 
lensing events, we can not accurately determine the 
amount of unblended lensed light because of blend- 
ing. In order to use these events when determining 
the microlensing optical depth, one approach could be 
to use the same technique as outlined above for pixel 
lensing. This would help to constrain a large number 
of classical events. Indeed this approach is desirable 
when we consider the effect of blending on measured 
optical depth. 

Aside from the aspects mentioned above, we note 
that, Griest and Hu (1992) and Dominik (1998) found 
that it was possible that galactic binary sources are 
fitted well with a single source with blending. As DIA 
does not suffer from blending this caveat is removed 
and the light curves should thus readily give unbiased 
results for binary sources. 

Again we remind the reader to the results of this 
analysis can be found in paper 2. 

7. SUMMARY 

We have outlined a detailed approach to the DIA 
scheme. A method of using multiple PSFs to ad- 
dresses the shortfalls of matching images with a single 
PSF was discussed. We noted how gradients in trans- 
mission and sky level across an image can be impor- 
tant when producing difference images. A solution 



using a spatially dependent photometric offsets and 
scaling was outlined. A method for reducing differ- 
ence images noise by combining observations to form 
the reference image was presented. 

The importance of using an accurate image reg- 
istration was stressed. The effect of differential re- 
fraction on difference images has been examined and 
is shown to be crucial. A method for compensating 
for differential refraction by offsetting stars relative 
to their colours was presented. 

A new method for determining the distribution of 
event time scales was discussed. This method would 
use existing HST colour-magnitude diagrams to de- 
termine the distribution of possible event sources with 
colour. The colours of events determined with DIA 
would then be used to determine a distribution of 
possible t values for each event. These would then be 
combined to determine the overall distribution of t 
values. From this lens mass distribution can be found 
in the usual way. 

We have demonstrated how the DIA photometry 
can make a large improvement in the quality of light 
curves over PSF photometry because of blending. We 
discussed how there are fundamental differences be- 
tween results from PSF and DIA photometry. Dif- 
ference images naturally provide unblended colours 
and centroid positions and lightcurves for microlens- 
ing events, while PSF results are usually blended to 
some extent. 
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Fig. 1. — Gaussian modelled profiles of three stars 
and the residuals when the two dashed profiles are 
differenced against the solid one. The short dashed 
lines corresponds to a centroid offset of 0.3". The 
long dashed lines corresponds to an offset of 0.07", 
the average offset of our data. 



Fig. 2. — Predicted differential refraction offsets be- 
tween stellar centroids for a range of star colours and 
airmasscs. Offsets arc simulated for three airmasses 
relative to a reference airmass of 1.01 and V-R colour 
of zero. The three airmass values are 1.9 (squares), 
1.5 (circles) and 1.2 (triangles). Results are for the 
MACHO Blue {B m ) passband. 
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Quotent Image pixel value 

Fig. 3. — The relationship between pixel values in the 
red divided by blue quotient image (colour map) and 
the standard MACHO passband colours B m and R m . 
Values are used to determine and correct for offsets 
cause by differential refraction for each pixel(star) in 
the image. 
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Fig. 4. — Difference in the offset induced by differ- 
ential refraction as a function of airmass, for a 5K 
(model) star and a 10K (model) star. The two MA- 
CHO Camera passbands B m and R m are shown. Off- 
sets in the blue passband are much bigger than the 
red because of the strong wavelength dependence of 
refraction. 



Fig. 5. — The effect of differential refraction on differ- 
ence images. Left image is a difference image without 
applying any correction for differential refraction ef- 
fects. The right image is the same image with our 
correction technique applied. The scale of the noise 
structures is much reduced although not completely 
removed. The two images are 100" x 100". The resid- 
ual object in the centre is due to a variable star 



Fig. 6. — Blending. Left: the position of an mi- 
crolensing event (arrowed) in an 11" x 11" ground 
based image Macho database image. Right: a drizzled 
HST image of the same region with reference stars (A 
- E) to approximately the same scale. Note the num- 
ber of objects blended within the ground based im- 
ages seeing disk. This result comes from an event in 
the LMC where the crowding is similar to the images 
analysed in this analysis. 




Fig. 7. — This figure shows the light curves for mi- 
crolensing event 96-BLG-16. The top two panels are 
the standard macho blue (B m ) and red (R m ) PSF 
photometry and microlensing fit residuals. The bot- 
tom two panels are the DIA photometry and resid- 
uals for the same initial data. The fits vary slightly 
because the PSF photometry is blended and the dif- 
ference image ph otometry is not. Further res ults are 
given in paper 2 ( Alcock ct al. (1999 papcr2)| ). 



Fig. 8. — Differences between microlensing light 
curves for a range of source fluxes. In the top panel 
we show three theoretical microlensing light curves 
for three values of source flux. In the lower panel we 
show the difference between the dashed curves and the 
solid curve. We can readily see that, to determine the 
differences between these light curves and hence true 
source flux, we require very small uncertainties and 
good sampling. The parameters are given in table [l| 



Table 1 

Parameters of Model Events. 



fs 


A 




i 


^max 


2.0 


2.40 


0.448 


85.4 


200 


1.5 


2.87 


0.365 


98.0 


200 


1.0 


3.83 


0.268 


120.8 


200 



Note. — Microlensing parameters 
for simulated events presented in fig- 
ure n 
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